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We present a measurement of moments of the inclusive hadronic-mass and -energy spectrum in 
semileptonic B — » X c l~V decays. This study is based on a sample of 232 million T(45) — > BB 
decays recorded by the BABAR detector at the PEP-II e + e~-storage rings. We reconstruct the 
semileptonic decay by identifying a lepton in events tagged by a fully reconstructed hadronic decay 
of the second B meson. We report preliminary results for the moments (m x ) with k — 1, ... ,6 and 
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(tlx) with k = 2, 4, 6 and = m\c A — 2AEx + A 2 , with mx the mass of the hadronic system, Ex 
its energy, and A a constant of 0.65 GeV, for different minimal lepton momenta between 0.8 and 
1.9GeV/e measured in the B- meson rest frame. These are predicted in the framework of a Heavy 
Quark Expansion (HQE), which allows the extraction of the total semileptonic branching fraction, 
the CKM-matrix element \V c b\, and the quark masses mi, and m c , together with the dominant 
non-perturbative HQE parameters. We find as preliminary results \V c t\ = (41.88 ± 0.81) • 10 -3 and 
m b = (4.552 ± 0.055) GeV/c 2 . 

Submitted to the 2007 Europhysics Conference on High Energy Physics, Manchester, England. 
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I. INTRODUCTION 

Measurement of moments of the hadronic-mass [1-5] 
and lepton-energy [4, 6, 7] spectra in inclusive semilep- 
tonic decays B X C £~V have been used to determine 
the non-perturbative QCD parameters describing these 
decays and the CKM matrix element \V c b\- 

Combined fits to these moments and moments of the 
photon-energy spectrum in B — > X s j decays [8-11] in 
the context of Heavy Quark Expansions (HQE) of QCD 
have resulted in precise determinations of \V c b\ and mb, 
the mass of the b quark. Specifically, they are reported to 
be \V cb \ = (42.0±0.2±0.7)-lCT 3 and m b = (4.59±0.03± 
0.03) GcV/c 2 in the kinetic mass scheme [12] and \V cb \ = 
(41.4±0.6±0.1)-10~ 3 andm 6 = (4.68 ± 0.03) GcV/c 2 in 
the IS scheme [13]. 

Lepton-energy moments are known with very good 
accuracy, but the precision of the hadronic-mass and 
photon-energy moments is limited by statistics. There- 
fore, we present here an updated measurement of the 
hadronic-mass moments (m x ) with k — 1, ... ,6 based on 
a larger dataset than previously used [2] . In addition we 
present measurements of the mixed hadron mass-energy 
moments (n|) with k = 2,4, 6 as proposed by Gam- 
bino and Uraltsev [14]. All moments are presented for 
different cuts on the minimum energy of the charged lep- 
ton. The mixed moments use the mass mx and the en- 
ergy Ex of the X c system in the B meson rest frame of 
B — > X c l~v decays, 

n 2 x = m 2 x c 4 - 2AE X + i 2 , (1) 

with a constant A, here fixed to be 0.65 GeV as pro- 
posed in [14]. They allow a more reliable extraction of 
the higher-order non-perturbative HQE parameters and 
thus they are expected to increase the precision on the 
extraction of \V cb \ and the quark masses nib and m c . 
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We perform a combined fit to the hadronic mass mo- 
ments, measured moments of the lepton-energy spec- 
trum, and moments of the photon energy spectrum in de- 
cays B — > X s -f. The fit extracts values for | V cb \ , the quark 
masses mj and m c , the total semileptonic branching frac- 
tion B(B — > X C £~T7), and the dominant non-perturbative 
HQE parameters. These are a\ and /ig, parameterizing 
effects at C(l/m 2 ), and p^ D and p\ Sl parameterizing ef- 
fects at 0(1/ ml). 



II. BABAR DETECTOR AND DATASET 

The analysis is based on data collected with the BABAR 
experiment [15] at the PEP-II asymmetric-energy e + e~ 
storage rings [16] at the Stanford Linear Accelerator Cen- 
ter between October 1999 and July 2004. 

The BABAR tracking system used for charged particle 
and vertex reconstruction has two main components: a 
silicon vertex tracker (SVT) and a drift chamber (DCH) , 
both operating within a 1.5-T magnetic field of a super- 
conducting solenoid. The transverse momentum resolu- 
tion is 0.47% at lGeV/c. Photons are identified in an 
electromagnetic calorimeter (EMC) surrounding a detec- 
tor of internally reflected Chercnkov light (DIRC), which 
associates Chercnkov photons with tracks for particle 
identification (PID). The energy of photons is measured 
with a resolution of 3 % at 1 GeV. Muon candidates are 
identified with the use of the instrumented flux return 
(IFR) of the solenoid. The detector covers the polar an- 
gle of 30° < < 140° in the center of mass (cm.) frame. 

The data sample consists of about 210.4 fb _1 , corre- 
sponding to (232 ± 3) x 10 6 decays of T (45) BB. We 
use Monte Carlo (MC) simulated events to determine 
background distributions and to correct for detector ac- 
ceptance and resolution effects. The simulation of the 
BABAR detector is realized with GEANT4 [17]. Simulated 
B meson decays are generated using EvtGen [18]. Final 
state radiation is modeled with PHOTOS [19]. 

The simulations of B — ► X c t~v decays use a pa- 
rameterization of form factors for B — > D*l~v [20], 
and models for B -> D£-V,D**£-V [21] and B -> 
DtxI-V,D*-eI-V [22]. 



5 



III. RECONSTRUCTION OF SEMILEPTONIC 
DECAYS 

A. Selection of Hadronic B-Meson Decays 

The analysis uses T(4S) — » BB events in which one 
of the B mesons decays to hadrons and is fully recon- 
structed (B reco ) and the semilcptonic decay of the re- 
coiling B meson (B rcco ii) is identified by the presence of 
an electron or muon. While this approach results in a 
low overall event selection efficiency of only a few per 
mille, it allows for the determination of the momentum, 
charge, and flavor of the B mesons. To obtain a large 
sample of B mesons, many exclusive hadronic decays are 
reconstructed [23]. The kinematic consistency of these 
-B r eco candidates is checked with two variables, the bcam- 
cnergy-substituted mass tties = \J s l^~ Vb anc ^ the en ~ 
ergy difference AE = Eb — y/s/2. Here y/s is the total 
energy in the cm. frame, ps and Eb denote the cm. 
momentum and cm. energy of the B roco candidate. We 
require AE = within three standard deviations, which 
range between 10 and 30 McV depending on the number 
of hadrons used for the reconstruction. For a given B reco 
decay mode, the purity is estimated as the signal frac- 
tion in events with toes > 5.27GeV/c 2 . For events with 
one high-momentum lepton with p| > 0.8GeV/c in the 
B-meson rest frame, the purity is approximately 78%. 



B. Selection of Semileptonic Decays 

Semilcptonic decays are identified by the presence of 
one and only one electron or muon above a minimum mo- 
mentum p\ min measured in the rest frame of the B reco ii. 
Electrons are selected with 94% average efficiency and a 
hadron misidentification rate in the order of 0.1%. Muons 
are identified with an efficiency ranging between 60% for 
momenta p = 1 GeV/c in the laboratory frame and 75% 
for momenta p > 2 GeV/c and a hadron misidentifica- 
tion rate between 1% for kaons and protons and 3% for 
pions. Efficiencies and misidentification rates are esti- 
mated from selected samples of electrons, muons, pions, 
and kaons. We impose the condition QbQt < 0, where Qi 
is the charge of the lepton and Qfc is the charge of the b- 
quark of the -B r cco ■ This condition is fulfilled for primary 
leptons, except for B°B° events in which flavor mixing 
has occurred. We require the total observed charge of 
the event to be |Q tat | = |Qb ioco + QB rccoil | < 1, allow- 
ing for a charge imbalance in events with low momentum 
tracks or photon conversions. In cases where only one 
charged track is present in the reconstructed X c system, 
the total charge in the event is required to be equal to 
zero. 



C. Reconstruction of the Hadronic System 

The hadronic system X c in the decay B — > X C £~T7 is 
reconstructed from charged tracks and energy depositions 
in the calorimeter that are not associated with the B reco 
or the charged lepton. Procedures are implemented to 
eliminate fake tracks, low-energy beam-generated pho- 
tons, and energy depositions in the calorimeter orig- 
inating from hadronic showers faking the presence of 
additional particles. Each track is assigned a specific 
particle type, either ^ , K^, or ir^, based on com- 
bined information from the different BABAR subdetectors. 
The four-momentum of the reconstructed hadronic sys- 
tem Px c is calculated from the four-momenta of the re- 
constructed tracks Pijrk, reconstructed using the mass 
of the identified particle type, and photons Pi :1 by 

Px c = Yh=i p i,trk +J2i=i p i,i- The hadronic mass 
rax is calculated from the reconstructed four-momenta 

as mx = \J P Xc - 

The four-momentum of the unmeasured neutrino P v 
is estimated from the missing four-momentum P m i SS = 
Pr(4S) ~ -Psreoo — Px c — Pi- Here, all four-momenta are 
measured in the laboratory frame. To ensure a well re- 
constructed hadronic system, we impose criteria on the 
missing energy, £?rniss > 0.5 GeV, the missing momentum, 
|Pmiss| > 0.5 GeV/c, and the difference of both quantities, 
|B miss — c|p m i SS || < 0.5 GeV. After having selected a well 
reconstructed -B reC o and having imposed the selection cri- 
teria on .Emiss — c|p miss |, 4.7% of signal decays and 0.3% 
of background decays are retained. 

Starting from a kincmatically well defined initial state 
additional knowledge of the kinematics of the semilcp- 
tonic final state is used in a kinematic fit to improve the 
overall resolution and reduce the bias of the measured 
values. The fit imposes four- momentum conservation, 
the equality of the masses of the two B mesons, and 
constrains the mass of the neutrino, P 2 = 0. The result- 
ing average resolutions in mx and n x are 0.355 GcV/c 2 
and 1.31 GeV 2 , respectively. The overall biases of the 
kinematically fitted hadronic system are found to be 
-0.096 GcV/c 2 and -0.08 GeV 2 , respectively We require 
the fit to converge, thus ensuring that the constraints are 
fulfilled. 

The background is composed of e + e~ — > qq(q = 
u, d, s, c) events (continuum background) and decays 
T(45) -> B+B~ or B°B° in which the B rcco candidate 
is mistakenly reconstructed from particles coming from 
both B mesons in the event (combinatorial background). 
Missing tracks and photons in the reconstructed hadronic 
system are not considered an additional source of back- 
ground since they only affect its resolution. The effect of 
missing particles in the reconstruction is taken care of by 
further correction procedures. To quantify the amount 
of background in the toes signal region we perform a 
fit to the toes distribution of the B [oco candidates. We 
parametrize the background using an empirical threshold 
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FIG. 1: tties spectrum of _B rC co decays accompanied by a 
lepton with p\ > 0.8GeV/c. The signal (solid line) and back- 
ground (red dashed line) components of the fit are overlaid. 
The crossed area shows the background under the -B rcC o sig- 
nal. The background control region in the rriBS sideband is 
indicated by the hatched area. 



function [24], 



diV 
dm ES 



oc m ES 



-x(i-* 2 ) 



(2) 



where x = m ES /m ES , max , m ES , max = 5.289 GeV/c 2 is the 
kinematic endpoint approximated by the mean cm. en- 
ergy, and x is a free parameter defining the curvature of 
the function. The signal is parameterized with a modi- 
fied Gaussian function [25] peaked at the i?-meson mass 
and corrected for radiation losses. The fit is performed 
separately for several bins in nix and n 2 x to account for 
changing background contributions in different m x or 
n x regions, respectively. The background shape is de- 
termined in a signal-free region of the to E s sideband, 
5.21 < m E s < 5.255 GeV/c 2 . Figure 1 shows the m E s 
distribution for p| > 0.8GcV/c together with the fitted 
signal and background contributions. 

Residual background is estimated from MC simula- 
tions. It is composed of charmless semileptonic decays 
B — > X U £~V, hadrons misidentified as leptons, secondary 
leptons from semileptonic decays of D^*\ mesons or 
r either in B°B° mixed events or produced in b — > ccs 
transitions, as well as leptons from decays of J/ip , and 
ip(2S). The simulated background spectra are normal- 
ized to the number of B reco events in data. We verify the 
normalization using an independent data control sample 
with inverted lepton charge correlation, QbQi > 0. 



IV. HADRONIC MASS MOMENTS 

We present measurements of the moments (m x ), with 
k = 1, ... 6, of the hadronic mass distribution in semilep- 
tonic B meson decays B — > X C £~V. The moments are 
measured as functions of the lower limit on the lepton 
momentum, p\ min , between 0.8 GeV/c and 1.9 GeV/c cal- 
culated in the rest frame of the B meson. 
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FIG. 2: Kinematically fitted hadronic mass spectra for mini- 
mal lepton momenta p\ > 0.8 GeV/c (top) and pi > 1.9 GeV/ c 
(bottom) together with distributions of combinatorial back- 
ground and background from non-BB decays (hatched area) 
as well as residual background (crossed area) 



A. Selected Event Sample 

The selected event sample contains about 21.5% back- 
ground. For p* > 0.8 GeV/c we find a total of 15085 ± 146 
signal events above a combinatorial and continuum back- 
ground of 2429 ± 43 events and residual background of 
1696 ± 19 events. For p* e > 1.9 GeV/c we find 2006 ± 53 
signal events above a background constituted of 271 ± 17 
and 248 ± 7 combinatorial and residual events, respec- 
tively. Figure 2 shows the kinematically fitted nix dis- 
tributions together with the extracted background shapes 
for p* > 0.8 GeV/c and p} > 1.9 GeV/c. 



B. Extraction of Moments 

To extract unbiased moments (m x ) , additional correc- 
tions have to be applied to correct for remaining effects 
that can distort the measured nix distribution. Con- 
tributing effects are the limited acceptance and resolution 
of the BABfiR detector resulting in unmeasured particles 
and in misreconstructed energies and momenta of par- 
ticles. Additionally measured particles not originating 
from the hadronic system and final state radiation of lep- 
tons contribute, too. We correct the kinematically fitted 
m x by applying correction factors on an event-by-event 
basis using the observed linear relationship between the 
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moments of the measured mass (m k x reco ) and moments 

of the true underlying mass (m| irae ). Correction func- 
tions are constructed from MC simulations by calculating 
moments (m x reco ) and (m x true ) in several bins of the 
true mass mx.true and fitting the observed dependence 
with a linear function. 

Studies show that the bias of the measured (m x reco ) 
is not constant over the whole phase space but depends 
on the resolution and total multiplicity of the recon- 
structed hadronic system, Nx c - Therefore, correction 
functions are derived in three bins of Nx c , three bins of 
-Bmiss — c|pmiss|, as well as in twelve bins of each with a 
width of 100 MeV/c. Due to limited number of generated 
MC events, the binning in Nx a and -Emiss — c|p m i ss | is 
abandoned for p| min > 1.7GeV/c. Overall we construct 
75 calibration functions for each order of moments. Fig- 
ure 3 shows examples of correction functions for the mo- 
ment {m x ) in three bins of p\ as well as in nine bins of 
-Emiss - c|p miss | and N Xa - 

For each event i the corrected mass m x cnIil) i is calcu- 
lated by inverting the linear function, 

fe _ ™X,reco.i ~ A i E nAss " c|p miss | , N Xc , fc, P* e ) 

m x>ca u b ,i - B{Emiss _ c \^ sa lN Xcl k,pl) 

(3) 

with A the offset and B the slope of the correction func- 
tion. Background contributions are subtracted by apply- 
ing weight factors Wi dependent on mx,reco to each cor- 
rected hadronic mass, whereby each weight corresponds 
to the fraction of signal events expected in the corre- 
sponding part of the mx,reco spectrum. This leads to 
the following expression used for the calculation of the 
moments: 

N ev 

S w i m X,calib,i 
( m x) = ' =1 Nev X Ccalib X C true . (4) 

i 

The factors C ca ub and Ctrue arc dependent on the order 
fc and minimal lepton momentum p\ min of the measured 
moment. They are determined in MC simulations and 
correct for small biases observed after the calibration. 
The factors C ca ub account for the bias of the applied cor- 
rection method and range between 0.985 and 0.996. For 
(m x ) we observe larger biases ranging between 0.897 
and 0.970 for the lowest p| min between 0.8GeV/c and 
1.2GcV/c, respectively. The residual bias correction fac- 
tor Ctme accounts for differences in selection efficiencies 
for different hadronic final states and QED radiation in 
the final state that is included in the measured hadron 
mass and distorts the measurement of the lepton's mo- 
mentum. The effect of radiative photons is estimated by 
employing PHOTOS. Our correction procedure results in 
moments which are free of photon radiation. The resid- 
ual bias correction Ctrue is estimated in MC simulations 
and typically ranges between 0.994 and 1.0f4. For the 
moments (m x ) and (m x ) slighly higher correction fac- 



tors are determined ranging between 0.994 and 1.023 as 
well as 0.994 and 1.036, respectively. 

This procedure is verified on a MC sample by applying 
the calibration to measured hadron masses of individual 
scmilcptonic decays, B — ► Dl~v, B — > D*i~v, four res- 
onant decays B — > D**£v, and two non-resonant decays 
B -> D^ixlv. Fi gure 4 shows the corrected moments 
{m x ) and (m x ) as functions of the true moments for 
minimal lepton momenta p\ > 0.8GcV/c. The dashed 
line corresponds to (m XtCaHb ) = (m x>true ). The cali- 
bration reproduces the true moments over the full mass 
range. 



C. Systematic Studies 

The principal systematic uncertainties are associated 
with the modeling of hadronic final states in scmilcp- 
tonic U-meson decays, the bias of the calibration method, 
the subtraction of residual background contributions, the 
modeling of track and photon selection efficiencies, the 
identification of particles, as well as the stability of the re- 
sults. The obtained results are summarized in Tables A.I 
and A. II for the measured moments {m x ) with fc = 1 . . . 6 
and selection criteria on the minimum lepton momentum 
ranging from p\ > 0.8 GeV/c to p\ > 1.9 GeV/c. 



1. Modeling of Signal Decays 

The uncertainty of the calibration method with re- 
spect to the chosen signal model is estimated by changing 
the composition of the simulated inclusive hadronic spec- 
trum. The dependence on the simulation of high mass 
hadronic final states is estimated by constructing correc- 
tion functions only from MC simulated hadronic events 
with hadronic masses mx,true,cut < 2.5GeV/c 2 , thereby 
removing the high mass tail of the simulated hadronic 
mass spectrum. The model dependence of the calibra- 
tion method is found to contribute only little to the total 
systematic uncertainty. We estimate the model depen- 
dence of the residual bias correction C trU e by changing the 
composition of the inclusive hadronic spectrum, thereby 
omitting one or more decay modes. We associate a sys- 
tematic uncertainty to the correction of the observed bias 
of the calibration method C ca ub of half the size of the ap- 
plied correction. 

We study the effect of differences between data and 
MC in the multiplicity and E mms — c|p m j SS | distributions 
on the calibration method by changing the binning of 
the correction functions. The observed variation of the 
results is found to be covered by the statistical uncertain- 
ties of the calibration functions. 
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FIG. 3: Examples of calibration fucntions for (nix) in bins of Nx c , E m i BB — c|p m i ss | and p|. Shown are the extracted (m X reco ) 
versus {m x true ) in bins of mx,true for 0.8 < p^ < 0.9GeV/c (•), 1.4 < p| < 1.5GeV/c (o), and p| > 1.9GeV/c (■). The results 
of fits of linear functions are overlaid as solid lines. A reference line with (rn X reco ) — (mx,true) i s superimposed (dashed line). 
There is only one calibration function with p} > 1.9GeV/c constructed but plotted for better comparableness in each bin. 



2. Background Subtraction 



Detector-Related Effects 



The branching fractions of background decays in the 
MC simulation are scaled to agree with current mea- 
surements [26]. The associated systematic uncertainty 
is estimated by varying these branching fractions within 
their uncertainties. At low Pg muV most of the studied 
background channels contribute to the systematic uncer- 
tainty, while at high p} min , the systematic uncertainty 
is dominated by background from decays B — > X u l~v. 
Contributions from Jji\) and if>(2S) decays are found to 
be negligible. 

The uncertainty in the combinatorial B rcco background 
subtraction is estimated by varying the lower and upper 
limits of the sideband region in the jties distribution up 
and down by 2.5MeV/c 2 . The observed effect is found to 
be negligible. 



We correct the MC simulation for differences to data in 
the selection efficiencies of charged tracks and photons, 
as well as identification efficiencies and misidentification 
rates of various particle types. The corrections are ex- 
tracted from data and MC control samples. 

The uncertainty of the photon selection efficiencies is 
found to be 1.8% per photon independent of energy, po- 
lar angle and multiplicity. The systematic uncertainty 
in track finding efficiencies is estimated to be 0.8% per 
track. We add in quadrature the statistical uncertainty of 
the control samples that depend on energy and polar an- 
gle of the track as well as the multiplicity of tracks in the 
reconstructed event. The misidentification of mesons 
as leptons is found to affect the overall normalization of 
the corresponding background spectra by 8%. 

While the latter two uncertainties give only small con- 
tributions to the total systematic uncertainty, the unccr- 



9 




- 






D % / 




,, * 




D /h 




& 


D* / 


M 


D n' 








w 




7i 





D. Results 



> [(GeV/c 2 ) 2 ] 



10 10 2 
<m*> [(GeV/c 2 ) 4 ] 



FIG. 4: Calibrated (•) and uncorrected (□) moments (m x ) 
(left) and (m x ) (right) of individual hadronic modes for min- 
imal lepton momenta p* t > 0.8GeV/c. A reference line with 
(■m x ,caiib) = (mx.true) is superimposed. 



The measured hadronic mass moments (m x ) with 
fe = 1 ... 6 as functions of the minimal lepton momen- 
tum pi min are depicted in Fig. 5. All measurements are 
correlated since they share subsets of selected events. Ta- 
bles A.I and A. II summarize the numerical results. The 
statistical uncertainty consists of contributions from the 
data statistics and the statistics of the MC simulation 
used for the construction of the correction functions, for 
the subtraction of residual background, and the deter- 
mination of the final bias correction. In most cases we 
find systematic uncertainties that exceed the statistical 
uncertainty by a factor of 1.5. 



MIXED HADRONIC MASS- 
ENERGY-MOMENTS 



AND 



tainty associated with the selection efficiency of photons 
is found to be the main source of systematic uncertain- 
ties. 



4- Stability of the Results 

The stability of the results is tested by dividing the 
data into several independent subsamplcs: B^ and B°, 
decays to electrons and muons, different run periods of 
roughly equal data-sample sizes, and two regions in the 

E m iss ^l^miss 

| spectrum, —0.5 < E miss — c|p miss | < OGcV 
and < E m i ss — c|p m iss| < 0.5 GeV, characterized by dif- 
ferent resolutions of the reconstructed hadronic system. 
No significant variations are observed. 

The stability of the result under variation of the se- 
lection criteria on -Emiss — c|p m iss| is tested by vary- 
ing the applied cut between |E m i ss — c|p m iss|| < 0.2 GeV 
and |-E m iss — c|p m i ss || < 1.4 GeV. For most of the mea- 
sured moments the observed variation is covered by other 
known systematic detector and MC simulation effects. In 
cases where the observed variation is not covered by those 
effects, we add an additional contribution to the system- 
atic uncertainty of the measurement that compensates 
the observed difference . 



5. Simulation of Radiation 

We check the impact of low energetic photons by re- 
moving EMC neutral energy deposits with energies below 
100 McV from the reconstructed hadronic system. The 
effect on the measured moments is found to be negligi- 
ble. 



The measurement of moments of the observable n 2 x , 
a combination of the mass and energy of the inclusive 
X c system, as defined in Eq. 1 , allow a more reliable 
extraction of the higher order HQE parameters /j, 2 , p? G , 
p 3 D , and p\g. Thus a smaller uncertainty on the standard 
model parameters \V c b\, in>bi an d m c could be achieved. 

We present measurements of the moments (n x ), (n x ), 
and (n x ) for different minimal lepton momenta be- 
tween 0.8 GeV/c and 1.9 GeV/c calculated in the B- 
meson rest frame. We calculate the central mo- 
ments ((n x — (n x )) 2 ), {{n 2 x — (n x )) 3 ), and the mo- 
ments ((n 2 x - 1.35 GeV 2 ) 2 ) and ((n| - 1.35 GeV 2 ) 3 ) as 
proposed in [14]. 

Due to the structure of the variable n 2 x as a differ- 
ence of two measured values, its measured resolution and 
bias are larger than for the mass moments and the sen- 
sitivity to -E m iss — c|p m iss| is increased wrt. to m x . The 
overall resolution of n 2 x after the kinematic fit for lep- 
ton momenta greater than 0.8 GeV/c is measured to be 
1.31 GeV 2 with a bias of -0.08 GeV 2 . We therefore intro- 
duce stronger requirements on the reconstruction quality 
of the event. We tighten the criteria on the neutrino ob- 
servables. The variable E m j ss — c|p m i ss | is required to 
be between and 0.3 GeV. Due to the stronger require- 
ments on -Emiss — c|Pmiss| the individual variables E m ; ss 
and [pmiss | have less influence on the resolution of the re- 
constrcutcd hadronic system. Therefore, the cuts on the 
missing energy and the missing momentum in the event 
are loosened to E m i ss > GeV and |p m i ss | > GeV/c, re- 
spectively, as they do not yield significant improvement 
on the resolution of n x , and do not increase the ratio of 
signal to background events. 

The final event sample contains about 22 % of back- 
ground events. The background is composed of 12 % con- 
tinuum and combinatorial background and 10 % decays 
of the signal B meson other than the semileptonic decay 
B — > X c i~V. Combinatorial and continuum background 
is subtracted using the sideband of the toes distribution, 
as described above. The residual background events, con- 
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FIG. 5: Measured hadronic mass moments (m\) with k = 1 . . . 6 for different selection criteria on the minimal lepton momen- 
tum Pi iln i n - The inner error bars correspond to the statistical uncertainties while the full error bars correspond to the total 
uncertainties. The moments are highly correlated. 



taining a correctly reconstructed i? rGC o meson, are sub- 
tracted using MC simulations. The dominant sources are 
pions misidcntificd as muons, B — > X U £~V decays, and 
secondary semileptonic decays of D and D s mesons. 

The measured n x spectra for cuts on the lepton mo- 
mentum at p* e > 0.8GcV/c and p\ > 1.9GeV/c 
are shown together with the backgound distributions in 
Fig. 6. We measure 7827 ± 108 (1278 ±42) signal events 
for p* > 0.8 (1.9)GeV/c, respectively. 



A. Extraction of Moments 

To extract unbiased moments (n x ), effects that dis- 
tort the n x distribution need to be corrected. These 
are the limited detector acceptance, resulting in a loss 
of particles, the resolution of measured charged particle 
momenta and energy depositions in the EMC, as well as 
the radiation of final-state photons. These photons are 
included in the measured X c system and thus lead to a 
modified energy and mass measurement of the inclusive 
system. In the case of radiation from the lepton, the lep- 
ton's measured momentum is also lowered w.r.t. its ini- 
tial momentum. The measured moments are corrected 
for the impact of these photons. 

As described before, we find linear relationships cor- 
recting the measured means (n Xrcco ) to the true means 
( n Jftme) described by first order polynomials. These 
functions vary with the measured lepton momentum, the 
measured E m { ss — c|p m i ss | , and the measured multiplic- 



ity of the inclusive X c system. The curves are therefore 
derived in three bins of -E^uss ~ c|p m i ss | and three bins of 
the multiplicity for each of the 12 lepton momentum bins 
of 100 MeV/c. We also find differences for events contain- 
ing an electron or a muon and therefore derive separate 
correction functions for these two classes of events. The 
measured n x value is corrected on an event-by-event ba- 
sis using the inverse of these functions: 

k _ n X,reco - M E miss - cjPmiss|, N Xc , k, p\) 

nX ' calib ~ B{E^- c\p^ ss \^x^Kp\) ■ ( > 

Here A and B are the offset and the slope of the calibra- 
tion function and differ for each order k = 2, 4, 6 and for 
each of the abovemcntioned bins. Figure 7 shows calibra- 
tion curves for the moments (n x ) (k = 2, 4, 6), integrated 
over all multiplicity bins and bins in E m i ss — c|p m i ss | , for 
three different bins of p* t . These calibration curves are 
extracted separately for events containing an electron or 
muon. Differences are mainly visible in the low momen- 
tum bin. 

To verify this calibration procedure, we extract the mo- 
ments of n x of individual exclusive B — > X C £~T7 modes 
on a MC sample and compare the calibrated moments 
to the true moments. The result of this study for the 
moments (n x ) is plotted in Fig. 8, confirming that the 
extraction method is able to reproduce the true moments. 
Small biases remaining after calibration are of the order 
of 1 % for (n 2 x ) and in the order of few percent for (n x ) 
and (n x ) and are corrected and treated in the systematic 
uncertainties. 
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FIG. 6: Spectra of n x after the kinematic fit together with 
distributions of combinatorial background and background 
from non-SB decays (hatched area) as well as residual back- 
ground (crossed area) for different minimal lepton momenta 
p* c > 0.8GeV/c (a) and p* c > 1.9GeV/c (b). 



Background contributions are subtracted applying n 2 x 
dependent weight factors Wi(n 2 x ) on an event-by-event 
basis, leading to the following expression for the mo- 
ments: 



JV e . 



E 



Xi 



'XI 



E Wi{n 2 x ) 

i=l 



xc(pt,k) 



(0) 



The bias correction factor C(pg,k) depends on the min- 
imal lepton momentum and the order of the extracted 
moments. It is derived on MC simulations and corrects 
for the small bias remaining after the calibration. 



B. Systematic Studies 



the neutral candidates and 0.8% for the charged tracks, 
corresponding to the systematic uncertainties of the ef- 
ficiency extraction methods. For the tracks we add in 
quadrature the statistical uncertainties from the control 
samples to the 0.8 % systematic uncertainty. The uncer- 
tainty arising from the differences between data and MC 
in the i? m i ss — c|p m i ss | distributions is evaluated by chang- 
ing the selected region of -Emiss — c|p m i ss | to [0.0,0.2] GeV 
and [0.0, 0.4] GeV. To evaluate the uncertainty due to 
the binning of the calibration curves in the multiplicity, 
we randomly increase the measured multiplicity used for 
the choice of the calibration curve by one with a probabil- 
ity of 5% corresponding to observed differences between 
data and MC. 

Smaller uncertainties arise from the unknown branch- 
ing fractions of the background decay modes. Their 
branching fractions are scaled to agree with recent mea- 
surements [26] and are varied within their uncertainties. 
The MC sample is corrected for differences in the iden- 
tification efficiencies between data and MC for various 
particle types. The uncertainty on the background due 
to pions misidentified as muons is evaluated by changing 
the MC corrections within the statistical uncertainties of 
these data control samples. While the background shape 
does not vary, the amount decreases up to 8 %. For the 
estimate of the uncertainty due to particle identification, 
we propagate this variation into the extracted moments. 

A similar variation procedure is applied for the branch- 
ing fractions of the exclusive signal modes, varying them 
several times randomly within 10% for the D*, 15% for 
the D, 50% for the individual D** modes and 75% for 
the non-resonant modes. The inclusive rate for the de- 
cays B — > X c l~v is conserved by rescaling all other 
modes. In addition, all D** (non-resonant) modes are 
scaled in common, again randomly within 50%, keeping 
the inclusive decay rate B — > X c l~v constant by rescal- 
ing the non-resonant (D**) modes only. Experimental 
uncertainties on the signal branching fractions are fully 
covered by these variations [26] . This dependence of the 
extraction method results in changes of the calibration 
curve and bias correction, however the impact on the 
moments measured on data is small. We conservatively 
add half of the bias correction remaining after calibration 
to the uncertainty related to the extraction method. 

The stability of the results has been tested by split- 
ting the data sample into several independent subsam- 
ples: B^ and B° , decays to electrons and muons, and 
different run periods of roughly equal data-sample sizes. 
No significant variations are observed. 



The main sources of systematic uncertainties have been 
identified as the simulation of the detector efficiency to 
detect neutral clusters. The corresponding effect from 
charged tracks is smaller but still contributes to the un- 
certainty on the moments. Their impact has been eval- 
uated by randomly excluding neutral or charged candi- 
dates from the X c system with a probability of 1.8 % for 



C. Results 

Figure 9 shows the results for the moments (n x ), 
(n x ), (n x ), and the central moments ({n x — {n 2 x )) k ) 
and ({n 2 x — 1.35 GeV 2 ) fc ) for k = 2 and 3 as a function 
of the p\ cut. The moments arc highly correlated due to 
the overlapping data samples. The full numerical results 
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FIG. 7: Examples of calibration curves for (fix) (k — 2,4,6) in bins of p\ , extracted separately for events B — > XeV (a)-(c) 
and B — > X/iV (d)-(f). Shown are the extracted (n| ireco } versus (rixtrue) m bins of n XtTue for 0.9 < p| < 1.0 GeV/ c (•), 
1.4 < p* t < 1.5GeV/c (□), and p\ > 1.9GeV/c (o) integrated over multiplicity and i? m i ss — c[p m i ss | bins. The results of fits of 
linear functions are overlaid as solid lines. Reference lines with {n Xrcco ) = (n^true) are superimposed (dashed lines). Please 
note the logarithmic scales in (b), (c), (e), and (f). 



and the statistical and the estimated systematic uncer- 
tainties can be found in Tables A. Ill - A. VII. A clear 
dependence on the lepton momentum selection criteria is 
observed for all moments, due to the varying contribu- 
tions from higher mass final states with decreasing lep- 
ton momentum. Statistical uncertainties on the moments 
(n k x ) arise from the limited data sample, the width of the 
measured distribution (n 2 ^) — (n^) 2 , and limited statis- 
tics on the MC samples used for the extraction of back- 
ground shapes, calibration curves, and bias correction. 
In most cases we obtain systematic uncertainties slightly 
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FIG. 8: Result of the calibration verification procedure for 
different minimal lepton momenta p* e > 0.8GeV/c (a) and 
Pi > 1.9GeV/c (b). Moments (n x ) of exclusive modes on 
simulated events before calibration (□) and after calibration 
(•) plotted against the true moments for each mode. The 
dotted line shows the fit result to the calibrated moments, 
the resulting parameters are shown. 



exceeding the statistical uncertainty. 



VI. DETERMINATION OF \V cb \ AND THE 
QUARK MASSES m b AND m c 



At the parton level, the weak decay rate for b 



cii 



can be calculated accurately; it is proportional to \V c b\ 
and depends on the quark masses, m& and m c . To relate 
measurements of the semileptonic B-meson decay rate to 
\V c b\, the parton- level calculations have to be corrected 
for effects of strong interactions. Heavy-Quark Expan- 
sions (HQEs) [27-29] have become a successful tool for 
calculating perturbative and non-perturbative QCD cor- 
rections [30-34] and for estimating their uncertainties. 

In the kinetic-mass scheme [14, 35-39], these expan- 
sions in 1/mb and a s (rrib) (the strong coupling constant) 
to order C(l/m^) contain six parameters: the running ki- 
netic masses of the 6— and c— quarks, TOb(/i) and m c (/i), 
and four non-perturbative parameters. The parameter 
fi denotes the Wilson normalization scale that separates 
effects from long- and short-distance dynamics. The cal- 
culations are performed for [i = 1 GeV [40] . We deter- 
mine these six parameters from a fit to the moments of 
the hadronic-mass and electron-energy [6] distributions 
in semileptonic B decays B — ► X c l~v and moments of 
the photon-energy spectrum in decays B — > X s ^ [10, 11]. 

In the kinetic-mass scheme the HQE to 0(l/mg) for 
the rate T$l of semileptonic decays B — > X c l~V can be 
expressed as [35] 
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FIG. 9: Measured moments (n 2 x ) (a), (n x ) (b), (n x ) (c), and the central moments {(n x — C) 2 ) with C = (n 2 x ) (•) and 
C = 1.35 GeV 2 (o) (d), and ((n 2 x - C) 3 ) with C = (n 2 x ) (•) and C = 1.35 GeV 2 (o) (e) for different cuts on the lepton 
momentum p| . The error bars indicate the statistical and the total errors, respectively. Please note the logarithmic scale on 
the j/-axis in plots (d) and (e). The moments are highly correlated. 
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The leading non-perturbative effects arise at 0(1 /to 2 ) 
and are parameterized by /i 2 (p) and [Tq (p) , the expecta- 
tion values of the kinetic and chromomagnctic dimension- 
five operators. At 0(1/to|), two additional parameters 
enter, p\)(p) and pf,s(/f), the expectation values of the 
Darwin and spin-orbit dimension-six operators, respec- 



tively. The ratio 



phase-space factor zq (r) = 1 — 8r + 8? 



I /ml enters in the tree level 



" 4 - 12r 2 lnr 



and in the function d(r) = 81nr + 34/3 - 32r/3 - 8r 2 + 
32r 3 /3 — 10r 4 /3. The factor 1 + A ew accounts for elec- 
troweak corrections. It is estimated to be 1 + A ew = 
(1 + a I V In Alz/mij) 2 = 1.014. The quantity A pert ac- 
counts for perturbative contributions and is estimated to 
be A pert (r,fx) « 0.908 [35]. 

The performed fit uses a linearized expression for the 
dependence of | V c b | on the values of heavy-quark param- 



eters, expanded around a priori estimates of these pa- 
rameters [35]: 



Wrl 



0.0417 



B 



civ 



1.55 



0.1032 t b 
x[l + 0.30(a s (m fc ) 

x [1 - 0.66(to 6 - 4.60) + 0.39(m c - 1.15) 



0.22)] 



0.013(^2 - 0.40) + 0.09O|, - 0.20) 
+0.05(/4 - 0.35) - 0.01(p| s + 0.15)]. (8) 

Here and m c arc in GeV/c 2 and all other parameters 
of the expansion are in GeV fc ; tb refers to the average 
lifetime of B mesons produced at the T(4S) and is given 
in ps. HQEs in terms of the same heavy-quark parame- 
ters are available for hadronic-mass, electron-energy, and 
photon-energy moments. Predictions for those moments 
are obtained from an analytical calculation. We use these 
calculations to determine |V^b|, the total semilcptonic 
branching fraction 23, the quark masses mj and to c , as 
well as the heavy-quark parameters p 2 ,, [Aq, p 3 Dl and p LS , 
from a simultaneous \ 2 fit to the measured moments and 
partial branching fractions, all as functions of the mini- 
mal lepton momentum p\ min and minimal photon ener- 
gies E 1<min . 
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A. Extraction Formalism 

The fit method designed to extract the HQE parame- 
ters from the moments measurements has been reported 
previously [12, 41]. It is based on a x 2 minimization, 



X 



= M, 



cxp 



Mi 



HQE 



M cxp — Mrqe 



(9) 



The vectors M exp and -Mhqe contain the measured mo- 
ments included in the fit and the corresponding mo- 
ments calculated by theory, respectively. Furthermore, 
the expression in Eq. 9 contains the total covariance ma- 
trix Ctot defined as the sum of the experimental, C cxp , 
and theoretical, Chqe, covariance matrices (see Section 
VIC). 

The total semileptonic branching fraction, B[B — > 
X c t~v\ is extracted in the fit by extrapolating mea- 
sured partial branching- fractions, B p * . (B — > X C £~T>), 
with p* ( > p\ min to the full lcpton energy spectrum. Us- 
ing HQE predictions of the relative decay fraction 



dTsLJ 771* 



Jo d M dE ? 



(10) 



the total branching fraction can be introduced as a free 
parameter in the fit. It is given by 



B(B -> XJ-V) 



B. 



a jB^XJ-V) 



R 



(11) 



The total branching fraction can be used together with 
the average B-meson lifetime tb to calculate the total 
semileptonic rate proportional to |V c b| 2 , 



SL 



B(B -> X c i~V) 



TB 



oc \V cb {< 



(12) 



By adding tb to the vectors of measured and predicted 
quantities, M exp and A/hqe, \V c b\ can be extracted from 
the fit as an additional free parameter using Eq. 8. 

The non-pcrturbative parameters fx G and p\ s have 
been estimated from B-B* mass splitting and heavy- 
quark sum rules to be p? G = (0.35 ± 0.07) GeV 2 and 
As = (-0.15 ± 0.10) GeV 3 [12], respectively Both pa- 
rameters are restricted in the fit by imposing Gaussian 
error constraints. 



B. Experimental Input 

The combined fit is performed on a subset of available 
moment measurements with correlations below 95% to 
ensure the invertibility of the covariance matrix. Since 
the omitted measurements are characterized by high cor- 
relations to other measurements considered in the fit they 
do not contribute significant additional information and 
the overall sensitivity of the results is not affected. All 
results are based on the following set of moment mea- 
surements, 27 in total: 



• Lcpton energy moments measured by BABAR [6]. 
We use the partial branching fraction B p * mea- 
sured for p\ > 0.6, 1.0, 1.5 GeV/ c and the moments 
(E e ) for p* > 0.6,0.8, 1.0, 1.2, 1.5GcV/c. The lep- 
ton energy moments (Ef) are used at the mini- 
mal lcpton momentum p\ > 0.6, 1.0, 1.5GeV/c and 
(Ef) at p* > 0.8,1.2GeV/c. 

• Hadronic mass moments are used as presented in 
this paper. We select the following subset for the 
fit: (m 2 x ) for p\ > 0.9, 1.1, 1.3, 1.5GeV/c and 



[in 



x) for P*e > °- 8 > 10, 1.2, 1.4GeV/c. 



• Photon energy moments measured in B — > X s j 
decays arc taken from [10] and [11]: (£7 ) for 
the minimal photon energy E 1 > 1.9, 2.0 GeV and 
(£ 2 ) for E 1 > 1.9 GeV . 

In addition we use tb = /o^o + (1 — fo)r± = (1.585 ± 
0.007) ps, taking into account the lifetimes [26] of neu- 
tral and charged B mesons, To and r±, and their relative 
production rates, /o = 0.491 ± 0.007 [26], the fraction of 
B°B° pairs. 



C. Theoretical Uncertainties 

As discussed in [12] and specified in [14] the following 
theoretical uncertainties are taken into account: 

The uncertainty related to the uncalculated pertur- 
bative corrections to the Wilson coefficients of non- 
perturbative operators are estimated by varying the cor- 
responding parameters /Lt 2 and n G by 20% and pf) and 
p\s by 30% around their expected values. 

Uncertainties for the perturbative corrections are es- 
timated by varying a s = 0.22 up and down by 0.1 for 
the hadronic mass moments and by 0.04 for the lepton 
energy moments around its nominal value. 

Uncertainties in the perturbative corrections of the 
quark masses mb and m c are addressed by varying both 
by 20MeV/c 2 up and down around their expected values. 

For the extracted value of \V c b\ an additional error of 
1.4% is added for the uncertainty in the expansion of the 
semileptonic rate Tsl [35, 39]. It accounts for remain- 
ing uncertainties in the perturbative corrections to the 
leading operator, uncalculated perturbative corrections 
to the chromomagnetic and Darwin operator, higher or- 
der power corrections, and possible non-perturbative ef- 
fects in the operators with charm fields. This uncertainty 
is not included in the theoretical covariance matrix Chqe 
but is listed separately as a theoretical uncertainty on 
\V cb \. 

For the predicted photon energy moments (E™), addi- 
tional uncertainties are taken into account. As outlined 
in [36], additional uncertainties of 30% of the applied 
bias correction to the photon-energy moments and half 
the difference in the moments derived from two differ- 
ent distribution-function ansatze have to be considered. 
Both contributions are added linearly 
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FIG. 10: The measured hadronic-mass and mixed moments (•/o), as a function of the minimal lepton momentum p^ min 
compared with the result of the simultaneous fit (solid line) and a previous measurement by the BABAR Collaboration (□) 
[2]. The solid data points mark the measurements included in the fit. The vertical bars indicate the experimental errors. The 
dashed lines correspond to the total fit uncertainty as obtained by converting the fit errors of each individual HQE parameter 
into an error for the individual moment. 



The theoretical covariance matrix Cr-qe is constructed 
by assuming fully correlated theoretical uncertainties for 
a given moment with different lepton momentum or pho- 
ton energy cutoff and assuming uncorrelatcd theoretical 
uncertainties for moments of different orders and types. 
The additonal uncertainties considered for the photon 
energy moments are assumed to be uncorrelatcd for dif- 
ferent moments and photon energy cutoffs. 

D. Results 

A comparison of the fit results for the hadronic-mass 
and mixed moments with the measured moments is 
shown in Fig. 10. The moments (m x ) and (m x } 
are not included in the fit and thus provide an unbiased 
comparison with the fitted HQE prediction. We find an 
overall good agreement, also indicated by \ 2 = 8 for 20 
degrees of freedom. The measured moments continue to 
decrease with increasing p\ min and extend beyond theo- 
retical predictions available for p\ min < 1.5 GeV/c. 

Comparing the measured moments (n x ) and 
{{n 2 x — (n x )) 2 } with predictions resulting from the 
presented fit, a good agreement is found. The calcu- 
lations used for the predictions of the mixed moments 
are currently missing p\ min -dependent perturbative 
corrections. The p* t min dependence of the perturbative 
corrections for those moments is however expected to be 
small [38]. 



The fit results for the standard model and HQE pa- 
rameters are summarized in Table I. We find as pre- 
liminary results \V cb \ = (41.88 ± 0.81) • 10~ 3 and m b = 
(4.552±0.055) GeV/c 2 . The results are in good agreement 
with earlier determinations [12, 13], showing slightly in- 
creased uncertainties due to the limited experimental in- 
put used in this fit. 

Figure 11 shows the A\ 2 = 1 contours in the (mi,, \ V c b\) 
and (mb,/i 2 ) planes. It compares the standard fit in- 
cluding photon energy moments, and a fit based on mo- 
ments from semileptonic B — > X c ^~77 decays only, clearly 
indicating the significance of the constraints from the 
B — > X s -f decays for both \V c b\ and m&. 



VII. SUMMARY 

We have reported preliminary results for the moments 
(m x ) with k = 1, . . . , 6 of the hadronic mass distribution 
in semileptonic i?-meson decays to final states containing 
a charm quark. In addition we have presented prelimi- 
nary results for a first measurement of the moments (n x ) 
for k = 2, 4, 6 with n x a combination of mass and energy 
of the hadronic system X c . The results for the mass mo- 
ments agree with the previous measurements [1-5] but 
tend in general to higher values, between 1% and 2% 
for (m x ) and (m%), respectively, relative to the previ- 
ous BABAR measurement [2]. The increased data sam- 
ple compared to the previous BABAR measurement led 
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TABLE I: Fit results with experimental and theoretical uncertainties. For |V c f,| we take an additional theoretical uncertainty 
of 1.4% from the uncertainty in the expansion of Tsl into account. Correlations coefficients for all parameters are summarized 
below the results. 





x 10 3 


m b [GeV/c^] 
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FIG. 11: A\ 2 = 1 contours for the fit results in the (mt, \V c b\) and (mt,/ij) planes comparing the results of the presented fit 
(black line) with those of a fit omitting the photon-energy moments (red dashed line). 



to significantly smaller statistical uncertainties which are 
smaller than the systematic uncertainties. 

We have made a combined fit in the kinetic scheme to 
the hadronic mass moments, measured moments of the 
lcpton-energy spectrum [6] , and moments of the photon 
energy spectrum in decays B X s -f [10, 11]. The com- 
bined fit yields preliminary results for \V c b\, the quark 
masses mf, and m c , the total semileptonic branching frac- 
tion B(B — > X C £~T>), and the dominant non-perturbativc 
HQE parameters in agreement with previous determi- 
nations. We obtain \V cb \ = (41.88 ± 0.81) • 10" 3 and 
m b = (4.552 ± 0.055) GeV/c 2 . 
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TABLE A.I: Results for the moments (m x ) with k = 1...3 for different cuts on the minimal lepton momentum p\ with 
absolute statistical and systematic uncertainties. Individual errors sources are specified due to modeling of the signal events, the 
calibration procedure, the background subtraction, detector efficiencies and resolution, and stability of moment measurements. 
Minimum lepton momenta cuts are given in GeV/c. Moments and uncertainties are given in (GeV/c 2 ) fe . 



k 


Pi ,min 


(m x ) 


O stat 


O sys 


Signal Model 


Ccalib 


BG subtr. 


Detector 


Stability 


1 


0.8 


2.0958 


±0.0083 


±0.0121 


0.0045 


0.0042 


0.0044 


0.0095 


0.0000 




0.9 


2.0920 


±0.0075 


±0.0107 


0.0039 


0.0040 


0.0042 


0.0082 


0.0000 




1.0 


2.0872 


±0.0072 


±0.0099 


0.0038 


0.0041 


0.0041 


0.0070 


0.0009 




f.f 


2.0796 


±0.0072 


±0.0093 


0.0036 


0.0035 


0.0041 


0.0066 


0.0000 




1.2 


2.0717 


±0.0075 


±0.0104 


0.0035 


0.0047 


0.0043 


0.0067 


0.0032 




1.3 


2.0661 


±0.0078 


±0.0128 


0.0032 


0.0054 


0.0045 


0.0067 


0.0077 




1.4 


2.0583 


±0.0081 


±0.0128 


0.0028 


0.0059 


0.0048 


0.0065 


0.0075 




1.5 


2.0518 


±0.0080 


±0.0121 


0.0025 


0.0063 


0.0053 


0.0071 


0.0045 




1.6 


2.0433 


±0.0089 


±0.0128 


0.0025 


0.0077 


0.0060 


0.0079 


0.0000 




1.7 


2.0378 


±0.0105 


±0.0162 


0.0024 


0.0075 


0.0073 


0.0080 


0.0091 




1.8 


2.0379 


±0.0139 


±0.0168 


0.0025 


0.0070 


0.0089 


0.0096 


0.0075 




1.9 


2.0350 


±0.0179 


±0.0225 


0.0020 


0.0098 


0.0121 


0.0121 


0.0107 


2 


0.8 


4.457 


±0.038 


±0.056 


0.022 


0.016 


0.018 


0.046 


0.000 




0.0 


4.430 


±0.032 


±0.048 


0.020 


0.014 


0.016 


0.038 


0.000 




1.0 


4.407 


±0.032 


±0.041 


0.019 


0.014 


0.015 


0.030 


0.000 




1.1 


4.368 


±0.031 


±0.039 


0.018 


0.011 


0.014 


0.029 


0.000 




1.2 


4.330 


±0.031 


±0.041 


0.017 


0.016 


0.014 


0.027 


0.015 




1.3 


4.299 


±0.032 


±0.051 


0.015 


0.018 


0.015 


0.027 


0.033 




1.4 


4.263 


±0.033 


±0.049 


0.014 


0.020 


0.016 


0.026 


0.031 




1.5 


4.231 


±0.031 


±0.045 


0.012 


0.021 


0.018 


0.027 


0.019 




1.6 


4.186 


±0.034 


±0.046 


0.012 


0.026 


0.020 


0.031 


0.000 




1.7 


4.157 


±0.040 


±0.056 


0.011 


0.024 


0.024 


0.031 


0.030 




1.8 


4.151 


±0.051 


±0.058 


0.011 


0.022 


0.029 


0.036 


0.025 




1.9 


4.128 


±0.065 


±0.077 


0.008 


0.035 


0.040 


0.045 


0.031 


3 


0.8 


9.67 


±0.15 


±0.21 


0.09 


0.05 


0.06 


0.17 


0.00 




0.9 


9.54 


±0.13 


±0.17 


0.08 


0.04 


0.05 


0.14 


0.00 




1.0 


9.45 


±0.11 


±0.14 


0.07 


0.04 


0.04 


0.10 


0.00 




1.1 


9.30 


±0.10 


±0.13 


0.07 


0.03 


0.04 


0.10 


0.00 




1.2 


9.16 


±0.10 


±0.13 


0.07 


0.04 


0.04 


0.09 


0.04 




1.3 


9.04 


±0.10 


±0.16 


0.06 


0.05 


0.04 


0.08 


0.11 




1.4 


8.92 


±0.10 


±0.15 


0.05 


0.05 


0.04 


0.08 


0.10 




1.5 


8.79 


±0.09 


±0.13 


0.05 


0.05 


0.05 


0.08 


0.06 




1.6 


8.62 


±0.10 


±0.13 


0.04 


0.06 


0.05 


0.09 


0.00 




1.7 


8.51 


±0.11 


±0.15 


0.04 


0.06 


0.06 


0.09 


0.07 




1.8 


8.47 


±0.14 


±0.15 


0.04 


0.05 


0.07 


0.10 


0.05 




1.0 


8.37 


±0.18 


±0.20 


0.03 


0.10 


0.10 


0.13 


0.06 
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TABLE A. II: Results for the moments (m\) with k = 4. . .6 for different cuts on the minimal lepton momentum p* t with 
absolute statistical and systematic uncertainties. Individual errors sources are specified due to modeling of the signal events, the 
calibration procedure, the background subtraction, detector efficiencies and resolution, and stability of moment measurements. 
Minimum lepton momenta cuts are given in GeV/c. Moments and uncertainties are given in (GeV/c 2 ) fe . 



k 


Pi .min 


(m k x ) 


0~stat 


V sys 


Signal Model 


Ccalib 


BG subtr. 


Detector 


Stability 


4 


0.8 


21.58 


±0.52 


±0.72 


0.30 


0.12 


0.21 


0.61 


0.00 




0.9 


20.98 


±0.41 


±0.58 


0.26 


0.09 


0.15 


0.48 


0.00 




1.0 


20.69 


±0.37 


±0.44 


0.26 


0.09 


0.12 


0.32 


0.00 




1.1 


20.17 


±0.33 


±0.40 


0.24 


0.06 


0.11 


0.29 


0.00 




1.2 


19.70 


±0.30 


±0.39 


0.23 


0.09 


0.10 


0.26 


0.13 




1.3 


19.28 


±0.29 


±0.46 


0.19 


0.10 


0.10 


0.23 


0.32 




1.4 


18.89 


±0.29 


±0.43 


0.17 


0.11 


0.10 


0.22 


0.29 




1.5 


18.45 


±0.27 


±0.36 


0.15 


0.12 


0.11 


0.23 


0.18 




1.6 


17.89 


±0.27 


±0.34 


0.13 


0.14 


0.12 


0.25 


0.00 




1.7 


17.51 


±0.31 


±0.37 


0.12 


0.13 


0.14 


0.25 


0.14 




1.8 


17.33 


±0.38 


±0.37 


0.11 


0.11 


0.17 


0.27 


0.09 




1.9 


16.98 


±0.46 


±0.49 


0.08 


0.24 


0.24 


0.34 


0.04 


5 


0.8 


49.83 


±1.97 


±2.47 


1.02 


0.40 


0.70 


2.10 


0.00 




0.9 


47.33 


±1.40 


±1.91 


0.90 


0.22 


0.47 


1.60 


0.00 




1.0 


46.40 


±1.23 


±1.39 


0.87 


0.21 


0.36 


1.00 


0.00 




1.1 


44.67 


±1.01 


±1.23 


0.80 


0.10 


0.29 


0.87 


0.00 




1.2 


43.22 


±0.93 


±1.15 


0.74 


0.19 


0.25 


0.75 


0.35 




1.3 


41.84 


±0.82 


±1.32 


0.62 


0.21 


0.24 


0.64 


0.92 




1.4 


40.69 


±0.82 


±1.21 


0.54 


0.23 


0.23 


0.59 


0.84 




1.5 


39.21 


±0.73 


±0.99 


0.45 


0.25 


0.25 


0.60 


0.53 




1.6 


37.46 


±0.73 


±0.88 


0.39 


0.30 


0.27 


0.68 


0.00 




1.7 


36.31 


±0.78 


±0.89 


0.35 


0.25 


0.32 


0.66 


0.26 




1.8 


35.68 


±0.96 


±0.87 


0.31 


0.20 


0.39 


0.68 


0.00 




1.9 


34.57 


±1.14 


±1.17 


0.23 


0.58 


0.53 


0.84 


0.00 


6 


0.8 


119.64 


±6.56 


±10.42 


4.19 


6.19 


2.28 


6.89 


0.00 




0.9 


109.40 


±4.49 


±7.80 


3.85 


4.13 


1.43 


5.19 


0.00 




1.0 


106.41 


±3.85 


±5.80 


3.83 


2.86 


1.03 


3.11 


0.00 




1.1 


100.70 


±3.22 


±4.73 


3.34 


1.87 


0.80 


2.66 


0.00 




1.2 


96.71 


±2.75 


±3.99 


2.95 


1.47 


0.60 


2.17 


0.00 




1.3 


92.56 


±2.45 


±3.77 


2.43 


1.07 


0.55 


1.67 


2.01 




1.4 


89.20 


±2.28 


±3.42 


2.10 


0.65 


0.54 


1.67 


1.94 




1.5 


84.54 


±1.99 


±3.03 


1.84 


0.58 


0.56 


1.63 


1.58 




1.6 


79.35 


±1.90 


±2.53 


1.62 


0.60 


0.61 


1.75 


0.00 




1.7 


76.01 


±2.01 


±2.58 


1.73 


0.46 


0.71 


1.68 


0.37 




1.8 


73.97 


±2.34 


±3.07 


2.41 


0.35 


0.87 


1.66 


0.00 




1.9 


70.67 


±2.71 


±4.76 


3.91 


1.34 


1.17 


2.04 


0.00 
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TABLE A. Ill: Results for {n x ) for k = 2,4,6 for all cuts pi . The systematic uncertainties are grouped in four categories 
having related sources: rec. efficiency is the sum of neutral and charged reconstruction efficiency differences data/MC, data/MC 
mismod. contains the errors from _E m i SB — c|p m i ss | differences and multiplicity differences, B bg. decays sums all contributions 
from the variation of the residual background component (including the fake lepton background), and signal model sums the 
impact of the variation of the signal decay branching fractions and the error related to the bias correction. 



k pi [GeV/c] 


(n L x) 




Cstat. 




<^"sys. 


rec. 
efficiency 


data/MC 
mismod. 1 


B 

)g. decays 


signal 
model 


2 0.8 


1.456 


± 


0.048 


± 


0.090 


0.054 


0.071 


0.010 


0.009 


0.9 


1.447 


± 


0.040 


± 


0.053 


0.038 


0.035 


0.010 


0.006 


1.0 


1.440 


± 


0.035 


± 


0.056 


0.042 


0.035 


0.008 


0.008 


1.1 


1.395 


± 


0.032 


± 


0.054 


0.038 


0.035 


0.006 


0.014 


1.2 


1.376 


± 


0.030 


± 


0.051 


0.035 


0.034 


0.004 


0.012 


1.3 


1.356 


± 


0.029 


± 


0.049 


0.032 


0.034 


0.004 


0.011 


1.4 


1.327 


± 


0.027 


± 


0.047 


0.032 


0.034 


0.004 


0.009 


1.5 


1.321 


± 


0.028 


± 


0.046 


0.030 


0.033 


0.004 


0.008 


1.6 


1.318 


± 


0.028 


± 


0.044 


0.027 


0.034 


0.004 


0.007 


1.7 


1.295 


± 


0.030 


± 


0.044 


0.027 


0.034 


0.005 


0.006 


1.8 


1.270 


± 


0.035 


± 


0.045 


0.028 


0.034 


0.006 


0.006 


1.9 


1.193 


± 


0.040 


± 


0.043 


0.021 


0.033 


0.007 


0.015 


4 0.8 


3.32 


± 


0.46 


± 


0.69 


0.34 


0.59 


0.08 


0.12 


0.9 


3.24 


± 


0.35 


± 


0.26 


0.17 


0.16 


0.07 


0.10 


1.0 


3.15 


± 


0.25 


± 


0.27 


0.19 


0.16 


0.05 


0.11 


1.1 


2.73 


± 


0.20 


± 


0.24 


0.13 


0.16 


0.04 


0.12 


1.2 


2.60 


± 


0.16 


± 


0.23 


0.11 


0.16 


0.01 


0.12 


1.3 


2.44 


± 


0.13 


± 


0.22 


0.09 


0.16 


0.02 


0.11 


1.4 


2.32 


± 


0.12 


± 


0.20 


0.07 


0.16 


0.01 


0.10 


1.5 


2.26 


± 


0.11 


± 


0.21 


0.08 


0.16 


0.01 


0.11 


1.6 


2.13 


± 


0.09 


± 


0.10 


0.06 


0.07 


0.00 


0.05 


1.7 


1.94 


± 


0.09 


± 


0.10 


0.05 


0.07 


0.00 


0.05 


1.8 


1.82 


± 


0.09 


± 


0.11 


0.06 


0.06 


0.00 


0.06 


1.9 


1.58 


± 


0.09 


± 


0.08 


0.05 


0.06 


0.00 


0.03 


6 0.8 


12.52 


± 


4.21 


± 


5.80 


2.83 


4.95 


0.63 


0.82 


0.9 


12.00 


± 


3.21 


± 


1.74 


1.13 


1.05 


0.55 


0.59 


1.0 


10.98 


± 


2.05 


± 


1.75 


1.22 


1.04 


0.35 


0.61 


1.1 


7.25 


± 


1.34 


± 


1.32 


0.63 


1.04 


0.22 


0.47 


1.2 


6.48 


± 


0.97 


± 


1.25 


0.54 


1.03 


0.09 


0.46 


1.3 


5.60 


± 


0.75 


± 


1.17 


0.35 


1.03 


0.08 


0.41 


1.4 


5.12 


± 


0.56 


± 


1.09 


0.17 


1.03 


0.03 


0.33 


1.5 


4.85 


± 


0.49 


± 


1.11 


0.19 


1.02 


0.02 


0.39 


1.6 


4.02 


± 


0.32 


± 


0.28 


0.16 


0.18 


0.02 


0.15 


1.7 


3.38 


± 


0.26 


± 


0.25 


0.12 


0.18 


0.01 


0.13 


1.8 


3.02 


± 


0.23 


± 


0.29 


0.14 


0.17 


0.01 


0.17 


1.9 


2.44 


± 


0.20 


± 


0.21 


0.11 


0.17 


0.02 


0.05 



TABLE A. IV: Results for ((n 2 x - (n 2 x }) 2 } for all measured cuts on p\ . 



Pi 


((n x - (n x )Y) 


Catat. 


C"stat.+sys. 


[GeV/c] 


[GeV 2 ] 






0.8 


1.20 


0.34 


0.57 


0.9 


1.15 


0.25 


0.28 


1.0 


1.08 


0.16 


0.21 


1.1 


0.78 


0.12 


0.17 


1.2 


0.71 


0.10 


0.16 


1.3 


0.60 


0.08 


0.14 


1.4 


0.56 


0.07 


0.13 


1.5 


0.52 


0.06 


0.14 


1.6 


0.39 


0.05 


0.08 


1.7 


0.27 


0.05 


0.08 


1.8 


0.21 


0.05 


0.09 


1.9 


0.16 


0.06 


0.09 



TABLE A.V: Results for ((n 2 x - 1.35 GeV 2 ) 2 } for all measured cuts on 



Pi 

[GeV/c] 


((n 2 x 


- 1.35 GeV 2 ) 2 ) 
[GeV 2 ] 


Cstat . 


(J stat.+sys. 


0.8 








0.40 


0.69 


0.9 




1. ID 




0.30 


0.35 


1.0 




1.09 




0.20 


0.28 


1.1 




0.78 




0.15 


0.22 


1.2 




0.71 




0.12 


0.20 


1.3 




0.60 




0.09 


0.18 


1.4 




0.56 




0.08 


0.16 


1.5 




0.52 




0.07 


0.16 


1.6 




0.39 




0.05 


0.06 


1.7 




0.27 




0.05 


0.04 


1.8 




0.21 




0.03 


0.05 


1.9 




0.19 




0.02 


0.02 



TABLE A. VI: Results for ((n| — {n 2 x )) 3 ) for all measured cuts on p| 

Pi ((n,x — Cstat. CTstat.+sys. 



[GeV/c] 


[GeV 3 ] 






0.8 


4.19 


2.38 


3.92 


0.!) 


3.99 


1.85 


1.99 


1.0 


3.33 


1.13 


1.36 


1.1 


1.26 


0.66 


0.84 


1.2 


0.96 


0.43 


0.67 


1.3 


0.68 


0.34 


0.59 


1.4 


0.57 


0.22 


0.53 


1.5 


0.51 


0.18 


0.53 


1.6 


0.18 


0.12 


0.17 


1.7 


0.18 


0.11 


0.16 


1.8 


0.19 


0.12 


0.19 


1.9 


0.16 


0.15 


0.21 



TABLE A. VII: Results for ((n% - 1.35 GeV 2 ) 3 ) for all measured cuts on 



Pi 


{(n x - 1.35 GeV 2 


) 3 ) fstat. 


(J stat.+sys. 


[GeV/c] 


[GeV 3 ] 






0.8 


4.57 


3.31 


5.58 


0.9 


4.33 


2.54 


2.85 


1.0 


3.62 


1.59 


2.04 


1.1 


1.36 


0.98 


1.33 


1.2 


1.01 


0.66 


1.08 


1.3 


0.69 


0.52 


0.94 


1.4 


0.53 


0.36 


0.82 


1.5 


0.46 


0.29 


0.80 


1.6 


0.14 


0.16 


0.22 


1.7 


0.13 


0.10 


0.16 


1.8 


0.14 


0.09 


0.16 


1.9 


0.08 


0.09 


0.15 



